GLOBAL BIFURCATION OF POSITIVE EQUILIBRIA IN NONLINEAR POPULATION 
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CHRISTOPH WALKER 

ABSTRACT. Existence of nontrivial nonnegative equilibrium solutions for age structured population models 
with nonlinear diffusion is investigated. Introducing a parameter measuring the intensity of the fertility, global 
bifurcation is shown of a branch of positive equilibrium solutions emanating from the trivial equilibrium. More- 
over, for the parameter-independent model we establish existence of positive equilibria by means of a fixed point 
theorem for conical shells. 



1. Introduction 

This paper deals with finding nonnegative equilibrium solutions to age and spatially structured population 
equations. In an abstract setting the problem reads: find a nontrivial function u : J — > Eq satisfying the 
nonlinear problem 

d a u + A(u, a) u + fj,(u, a) u = , a 6 J \ {0} , (1.1) 

u(0) = [ 0(u, a) u(a) da , (1.2) 
Jo 

where Eq is an ordered Banach space with positive cone Eq and J := [0, a m ) with a m £ (0, oo] denotes the 
maximal age. Equations ( 11.1b . (11.21 i arise naturally when considering equilibrium (i.e. time-independent) 
solutions to population models, where the function u represents the density of a population of individuals 
with age a £ J whose evolution is governed by death and birth processes according to the density dependent 
death modulus fj,(u, a) and birth modulus j3(u, a), respectively. The term A(u, a) is (for u and a fixed) a 
linear unbounded operator A(u, a) : E\ C Eq — > Eq defined on some common subspace E\ of Eq and 
models spatial movement of individuals. The full non-equilibrium equations involve an additional time 
derivative in (II. j} , (11.2b . Such age-structured equations have been studied since long ago (see [29 30 1 and 
the references therein), in particular in situations where spatial movement is neglected (i.e. A = 0) or when 
A does not depend on the density u itself (see fl7][T6][18 23 24] and the references therein). Less seems to 
be known about models involving nonlinear age-dependent diffusion (however, see I151 [151 [1711261 1271 ). 

To understand the asymptotic behavior of the time evolution of structured populations a precise knowl- 
edge about equilibrium solutions (i.e. solutions to ( II. lb . (11.2b ) is needed. In the present paper we focus on 
such solutions for the nonlinear age-dependent case A = A(u, a). Clearly, u = is a solution to ( II. lb . 
( II. 2b and thus the aim is to give conditions for existence of nontrivial solutions. Moreover, since u in dl.lb , 
( II. 2b represents a density, any solution should be nonnegative an thus, in the abstract setting, belong to the 
positive cone Eq. 

To shorten notation we introduce an operator A as 

A(u, a) := A(u, a) + /j,(u, a) . (1.3) 

Suppose that for u fixed, the map a A(u, a) generates on the Banach space Eq a parabolic evolution 
operator IT u (a, a), < a < a < a m . Then an easy - but fundamental - consequence of properties of 
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evolution operators is that any solution u to dl.lt , ( 11. 2t must satisfy the relation 

u(a) = U u (a, 0)u(0) , aeJ, u(0) = Q{u)u{0) , (1.4) 
where the linear operator Q(u) on Eq is (for u fixed) given by 

Q(u) :— / /3(u, a)U u (a, 0) da . 
Jo 

Roughly speaking, Q(ii) contains information about the spatial distribution of the average number of off- 
spring per individual over the entire lifespan of the individual. If spatial movement is neglected, that is, if 
A = and hence n tl (a, a) = e~ ft i- l ( u ^ r ) dr ) then Q(u) is simply the net reproduction rate (see [29|), and 
for any solution u to dl.lt , d 1 .2b with A = this number Q(u) necessarily equals 1 according to ( 11. 4t . If 
spatial movement is included, then (11.4b implies that ti(0) is (if nonzero) an eigenvector to the eigenvalue 1 
of the operator Q(u). 

In Section|2]we suggest a bifurcation approach to establish positive solutions emanating from the trivial 
solution u = 0. We introduce a bifurcation parameter n, which determines the intensity of the fertility 
without changing its structure, by setting 

nb(u,a) := f3(u,a) , (1-5) 
where b is normalized such that the spectral radius r(Qo) of the bounded linear operator 

Qo ■= f b(0,a)n (a,0)da 
Jo 

satisfies r(Qo) — 1. Hence r(Q(0)) = nr(Qo) = n and n thus represents the "inherent net reproduction 
rate at low densities" (technically when u = 0). If Qo is a compact positive operator, then 1 is an eigenvalue 
of Qo and there is a distribution for which spatial movement, birth, and death processes balance each other 
if these processes are governed according to A(0, •), 6(0, •), and /i(0, •). 

In ll28l a local bifurcation result was shown and a set of n-values around the critical value n = 1 was 
provided for which ( ll.lt , d 1 .2b subject to dl.5t possess nontrivial positive solutions. The aim of this paper 
is to extend these results to prove a global bifurcation phenomenon. More precisely, if the operator A(u, a) 
admits a suitable decomposition A(u,a) = Aq(cl) + A* (it, a) with A* being of "lower order" (see Re- 
mark !2.ll and d2.61 >) we use Rabinowitz's alternative [22 1 to show that there is an unbounded continuum of 
nontrivial positive solutions (ri, u) to ( ll.lt . (II .2b subject to d 1 .5b . Furthermore, we characterize the set of 
n-values more detailed in some cases and give an example demonstrating that the assumptions imposed are 
quite natural. We shall point out that our results and methods were inspired by [9 1, where global bifurcation 
for population models neglecting spatial structure from the outset were investigated. More results on bifur- 
cation for age structured equations with A = can be found in [8] [9] [101] and, respectively, in lfl2l [T3 1 for 
linear and age-independent A. We also refer to |fT9"l l20l |2T1 129 | and the references therein for equilibrium 
solutions for age structured equations in general. 

The subsequent section is then devoted to a different approach for establishing solutions to the parameter- 
independent problem (ll.lt . (II .2b not assuming a particular decomposition of A. This approach covers 
"fully" quasilinear problems and is more or less independent of the previous considerations. Note that the 
form of the solution in (11. 4t allows one to interpret (u, B) with B := u(0) as a fixed point of the map 

(u,B) ^ (ll u (;0)B,Q(u)B) . 

In Section [3] we extend an argument of ||20l Thm.l] (see also [29 Thm.4.1]) for non-diffusive population 
equations which is based on a fixed point theorem for conical shells [ 1 1. The use of such a theorem prevents 
hitting the trivial solution u = (together with B = being obviously a fixed point of the map above). 
We thus prove existence of nontrivial positive solutions for (ll.lt . (II .2b under fairly general assumptions. 
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Loosely speaking, nontrivial positive solutions exist provided that the spectral radius of Q(u) for small pop- 
ulations with density u satisfies r{Q(u)) — r(Q(Q)) = 1 and is an eigenvalue with a common eigenvector 
B for all Q(u), and provided that, in addition, for large populations densities there holds r(Q(u)) < 1. 
Thus relevant equilibrium solutions exist if small populations do not affect the spatial distribution of net 
reproduction rate and large populations have a spatial net reproduction rate not exceeding 1 . 

We conclude the introduction with some notation being used in the following. If E and F are Banach 
spaces we write C(E, F) for the set of linear bounded operators from E to F, and we put C(E) := C(E, E). 
The subset thereof consisting of compact operators is denoted by K.(E, F) and JC(E), respectively. We write 
r(A) for the spectral radius of an operator A G C(E). For an ordered Banach space E we let C+ (E) denote 
the positive linear operators and K.+ (E) is the set of compact positive linear operators. Next, £is(E, F) 

d 

stands for the set of topological isomorphisms E — > F. By E <—* F we mean that E is densely embedded 

in F and E CJ ^ F stands for a compact embedding of E in F. If E ^ F we let Ti(E, F) denote the set of 
all negative generators of strongly continuous analytic semigroups on F with domain E. Moreover, given 
lj > and k > 1 we write A G H(E, F; k, oj) if A G C(E, F) is such that lu + A e Cis(E, F) and 

~< nn X t A)U ! F n <"» ^A> W , ueE\{0}. 

K ~ |A| ||tt|| F + ||u||jS ~ 



Note that 



H(E,F)= |J H(E,F;k,u) 



K>1 
UJ>0 



We refer to [4| for more details. 

Given open subsets X C E, Y C F and another Banach space G we mean by / G C°' P (X x Y, G) for 
p > a continuous map / : X x Y — > G such that f(x, •) is p-H61der continuous for each x E X. Then 

r,, o \\f(x,Vi) - f(x,y 2 )\\ G 

P ' VUV2ZY \\yi-y2\\F 

We let Cb{E, F) denote the continuous functions from E to F being bounded on bounded sets. 



General assumptions. Throughout the paper we assume that Eq is a real Banach space ordered by a closed 

convex cone E+ and E x £ for some Banach space E\. We fix for each 8 G (0, 1) an admissible 

interpolation functor (•, that is, an interpolation functor (■, -)g such that E\ := (i?o, £7i )e» ■ Note 

that Eg CJ —* E$ for < -d < 9 < 1 (see J4] I.Thm.2.1 1.1]). The interpolation spaces Eg are given their 
natural order induced by the cone Eg := Eg n i?^. We fix a m G (0, 00] and set J := [0, a m ). Observe 
that a m = 00 is explicitly allowed. 

2. Global Bifurcation of Positive Equilibria 

We focus our attention on the parameter dependent problem ( 11.11 ). ( 11.2b subject to ( 11.51 ) in a more general 
framework. More precisely, we look for solutions (n, u) to problems of the form 

d a u + A(u, a) u = , a G J , 

u(0) = nf(u) . 

The linear unbounded operator A (it, a) : E\ C -Bo — » £0 (for it and a fixed) and the operator £ with 
i{u) G £q and £(0) = are supposed to satisfy some technical assumptions specified later on. The main 
example for the latter we have in mind is, of course, 

£(u) = / b(u,a)u(a) da . 
Jo 
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Clearly, the branch (n, u) = (n, 0), n £ R, consists of (trivial) solutions to ( 12. It . Our aim is to prove that 
another unbounded branch of nontrivial positive solutions (n, u) (i.e. u(a) G Eq for a £ J and it ^ 0) 
bifurcates from the trivial branch at some critical value, which we may assume to be n = 1 under a suit- 
able normalization. Imposing maximal i p -regularity for (a part of) the operator A we will show that the 
nontrivial branch (n,u) is unbounded in R+ x (L+(J,Ei) n Wp(J,E )) for some p G (1, oo) fixed. The 
result is inspired by the work of |9| for the non-diffusive case A = and is a consequence of Rabinowitz's 
alternative El . The application of this alternative requires some compactness in appropriate spaces of 
the operators involved, which is guaranteed, for example, by a generalized Aubin-Dubinskii lemma (see 
Remark |2~T| and Lemma|23Tl. 



To state the precise assumptions let us fix p £ (1, oo) and introduce the spaces 
E := L p {J, E ) , Ei := L p {J, E 1 ) n W p x {J, E ) . 

Recall the embedding 

Ei ^BUC{J,E q ) , (2.2) 

where BUC stands for bounded and uniformly continuous and E q := (Eo, Ei)^ tP with (•, p being the 
real interpolation functor for ? := := 1 — 1/p. The trace 'yu := u(0) is thus well-defined for u £ Ei. 
We further fix Banach spaces Fj., F2, F3, and F4 such that 

Ei^F, , j = 1,2, 3, 4, (2.3) 

and first remark the following: 

Remark 2.1. Given a £ [0, 1) and s £ [0, 1 — a), Sobolev spaces Wp(J, E a ) are appropriate choices for 
Fj to satisfy ( 12. 31 ). This follows from the compact embedding E\ E a and a generalized Aubin-Dubinskii 
lemma [5 , Thm. 1 . 1 ] ( together with a diagonal sequence argument if J = M. + ). 

For the nonlinear operator A in (12. U we then shall assume a decomposition of the form 

A(u,a) = Ao(o) + A*(w,a) , (2.4) 

where Ao is an age-dependent parabolic operator and the nonlinearity of A in it is contained in a "lower 
order perturbation" A*. To be more precise we suppose for the linear part Ao that 

Ao £ Loo(J, C(E%, Eo)) generates a positive parabolic evolution operator 

IIo(a, a), < a < a < a m , on Eq with regularity subspace E\ and (2.5) 

possesses maximal L p -regularity, that is, (d a + Ao, 7) £ £is(Ei, Eo x E q ) , 

while for the nonlinear part A* we assume that 

A, g C(Fi, C(¥ 2 , E )) with A*(0, •) = . (2.6) 

We also assume £ in (12. U admits a decomposition 

t(u) = £ (u) + l*(u) (2.7) 

with linear part 

£ a £ £(Ei, £0) + £(F 3 , £: f ) for some 1? G (?, 1] , (2.8) 

and nonlinear part 

4 6a(Ei,£,) + C(F 4l £ t ) with ||4(u)|U t =o(||u|| El ) as ||u|| El -> . (2.9) 
A consequence of ( 12.5b is that for each datum (it , /) G E^ x Eo the problem 

d a u + A (a)u — /(a) , a £ J , u(0) = u 
possesses a unique solution u £ Ei given by 

u(a) = IIo(a,0)u + (tfo/)(a), 
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satisfying for some number cq > independent of / and u° 

Nk<c&(ll«V + ll/lk>) » 

where 

(Kb/) (a) := / n Q (a, a)f(a) da , a G J , / G E . 
Jo 

Therefore, 

n (-,0) e£(£ s ,Ei) , A'cG^EcEi), (2.10) 

and (O, (O thus imply 

[/ h+ 4(ff /)] e£(E ,S ? ) . (2.11) 

We also note that 

Qo G K(E f ) for Q w := € (Ho(-,0H , w e E t , (2.12) 

which is a consequence of ( 12.101 1 and either E$ E^ if £ G £(Ei,.E^) or (12.3b if £o G £(F 3 , 
Without loss of generality we may assume that £q is normalized such that the spectral radius of Qo G IC(E^) 
equals 1, that is, r(Qo) = 1. We first consider the linearization of d2.lt around u = 0. 

Lemma 2.2. Suppose (I2.3b -( |2~9l ) one/ /ef r(Qo) = 1- Then, for each (hi, h 2 ) G -E ? x Eo, f/ie problem 

d a u + A (a) u = h 2 (a) , a G J , u(0) — 7^0 ( u ) = hi (2.13) 
admits a unique solution u = /la) G Ei given by u — Hq(- , 0)w(hi, h<i) + Koh 2 , where 

1 \-Vi 



w(h u h 2 ) := ^1- -Q j ^-£ {K k 2 ) + h 1 jeE i . (2.14) 
The solution operator S belongs to C(E^ x Eo, Ei). 

Proof. By the previous observations problem ( 12.131 ) is, for any given (hi, h 2 ) G E q x Eo, equivalent to 

u(a) = n (a,0)u(0) + (K h 2 ){a) , aeJ, 

u(0) = ^Q u(0) + ^£ {K h 2 ) + hi . 

Taking 2 > r(Q ) into account, the latter equality entails ( 12. 14b , and defining 

S(h u h 2 ) := n (-, 0)w(hi, ha) + K h 2 , 

we derive S G C(E, x E , Ei ) from J2TT0b and d2TTTT >. □ 

Lemma l2~2l allows a reformulation of problem ( 12. U in terms of the operator 5: writing n = A + 1/2, a 
function it G Ei solves ( 12.11 ) if and only if 

u = XLu + H(X,u) , (2.15) 

where 

Lu := S(£ (u),0) , H(X,u) := S((X + l/2)f,(u), -A*(u, •)«) ■ 

The maps L and enjoy the following properties. 

Lemma 2.3. Suppose (l2T3T>-d2T9T> a/id /ef r(Q ) = 1. ^/zen L G /C(Ei). Moreover, # G C(M x Ei,Ei) is 
compact and \\H (A, w)IIei = °(|| u l|iEi) a* |M|ei - * uniformly on bounded X intervals. 
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Proof. Clearly, Lemma O (12731 , and <|27J) ensure L E £(E X ). Using ^ £ ? if £ £ £(Ei,J3tf) 
or Ei F3 if Iq € £(F3, £^), the compactness of L is obvious. Next, Lemma [2721 together with (12.6b 
and ( 12.91 ) imply JI e C(l x Ei, Ei). As for its compactness we note that if (uf) is a bounded sequence 
in Ei, then (£*(uj)) is relatively compact in either because {£*{uj)) is bounded in E$ E s if 
4 G C b (Ei, £#) or because Ei F 4 if 4 G C(F 4 , £?J. Next, (IQT i and (|2~6l > ensure that (A*(uj)uj) is 
relatively compact in E . Hence the compactness of H follows from Lemma [2721 Finally, the last assertion 
is a consequence of ( 12.61 l. ( 12.9b , and again Lemma [2~2l □ 

Remark 2.4. Alternatively to d2.6b we could have assumed A* G C(Ei, £(Ei,Eo)) w/fn A*(0, •) = is 
smc/i that (it 1— > A*(it)it) G C(Ei, Eo) ts compact. 

To problem ( 12.15b we may now apply Rabinowitz's alternative l22l . Recall that the characteristic values 
of a linear operator are the reciprocals of its real nonzero eigenvalues. A continuum in R x Ei is a closed 
connected subset thereof, and it meets infinity if it is unbounded. 

Proposition 2.5. Suppose (12.3b - d2~9l i and let r(Q ) = 1 be a simple eigenvalue of Qq G IC(E^) with 
eigenvector B G E^. Then there exists a maximal continuum £ in R X Ei consisting of solutions {n, it) to 
(12. lb wtf/z it ^ z/n ^ 1 ant/ (1, 0) G £ -For (n, it) G £ near (1, 0) we nave 

it = eIIo(-, 0)B + u*(e) w/fn ||u*(e)|| El = o(e) as e^O. (2.16) 

Tne continuum £ satisfies the following alternative: either £ meets infinity or it meets a point (n, 0) with 
jj, = n — 1/2 being a characteristic value of L. 

Proof. First note that u = pLu with iigEj and /ieRis equivalent to 

u = n (-, 0)«(0) , u(0) = (/1 + 1/2)Q «(0) (2.17) 

owing to Lemma |275] and (12.2b . Thus, /x + 1/2 is a characteristic value of Qo E if an d on ly if 

/1 is a characteristic value of L G /C(Ei). Additionally assuming /x + 1/2 to be a simple characteristic 
value of Qq we claim that p is a simple characteristic value of L. For, let it G ker((/iX — l) 2 ) C Ei and 
define w := (/xL — l)it G ker(/xL — 1) C Ei. Then ( 12.17b ensures i> = IIo(-, 0)u(0) with v(0) belonging 
to ker(l — (/x + l/2)Qo). The characteristic value p + 1/2 of Qo being simple, we deduce v(0) — r£o for 
some r G R and £0 G E q with ker(l — (p + l/2)Qo) = span{£o}- Hence 

u = rn (-,0)£o, (2.18) 

and we aim for r = 0. Clearly, from (12.17b we have no(-,0)£o € ker(/iL — 1), whence from ( 12.18b 
u = L(pu — rfj,Ho(-, 0)£o)- Lemma [2721 then entails 

(1 - (p + l/2)Q )u(0) = -r/xQ £o = ft • 

Consequently, 

r/i^o G rg(l - O + 1/2)Q ) D ker(l - (p + 1/2)Q ) = {0} 

since /x + 1/2 is a simple characteristic value of the compact operator Qo- We conclude r — as desired 
because p = is impossible owing to the fact that 1/2 is no characteristic value of Q since r(Q ) = 1. 
But then v = by d2.18b and so ker((/xZ/ — l) 2 ) C ker(/xL — 1). Therefore, /i is indeed simple for L 
provided p, + 1/2 is simple for Qo- In particular, p = 1/2 is a simple characteristic value of L due to the 
assumption on r(Qo). 

Taking Lemma [2731 into account we may apply Theorem 1.3 from j22l to conclude the existence of a 
maximal continuum £ in 1 x Ei with (1,0) G £ such that (n, u) G £ solves it = XLu + H(X, it) with 
A = n — 1/2, where it ^ if n ^ 1, and £ either meets infinity or meets a point (n, 0) with a characteristic 
value /t = h — 1/2 of L different from 1/2. Finally, ll22l Lem.1.24] implies ( 12.16b and the statement 
follows. □ 
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Imposing further conditions on A and I we now prove that a global branch of positive solutions (n, u) to 
(12. U exists emanating from the trivial branch (n, 0), n G R at the critical point (1, 0). To prove this result 
we suppose that 

for each u G Ei, A(u, ■) generates a positive parabolic evolution operator 

(2.19) 

H u (a, a), < a < a < a m , on Eq with regularity subspace E\ . 
We also assume that 

there is with £*(u) — £*(u, u) and £«(0, •) = such 

that Q u G /C + (£ , ? ) for each u G Ei , where (2.20) 

q u w :=^ (n„(-,oH +f»(!i,n„(.,ow,w e £ ? . 

Note that this definition of Q u is consistent with ( 12.121 ). Let 'mt(E^) denote the interior of the positive cone 
E+. Then we assume further that 

for each u G Ei, any positive eigenvector to a positive eigenvalue of Q u belongs to mt(E^) . (2.21) 

This last assumption is crucial for positivity of solutions but not too restrictive in applications as noted in 
the following remark (see also Example 12. 10b . 

Remark 2.6. If mt(££") and Q u G JC+(E q ) is irreducible (e.g. if strongly positive), then the Krein- 
Rutman theorem 1111 Thm.12.3] ensures that the spectral radius r(Q u ) > is a simple eigenvalue of Q u 
with an eigenvector belonging to int(E^), and it is the only eigenvalue with positive eigenvector. Thus 
( 12. 2U holds in this case. 

Theorem 2.7. Suppose ( 12.31 l-( |2~9l ), and ( |2.19l )-( |2.21| i. Further let r(Qo) — I be a simple eigenvalue of 
Qo G /C+(£J<j) with eigenvector B G mt(E^) and suppose Qo has no other eigenvalue with positive 
eigenvector. Then there is a contiunuum C + in R + x Ej 1 " of positive solutions to d2.ll) connecting (1,0) with 
infinity. 

Proof. Let £ denote the maximal continuum of solutions to ( 12. U provided by Proposition [23] Clearly, if 
(n, u) G £, then u = II„(-, 0)u(0) by (O), j2l5] |. and ( |2~T9] l. Thus u G E+ provided u(0) G E+. Also 
note that (0, u) G £ would imply u = and is thus impossible. Due to B G 'mt(E^) it follows from d2.2l) 
and ( 12.161 ) that for (n, it) G £ near (1, 0) we have for sufficiently small e > 

-u(0) = B + -fuJe) G . 

e e 

Since /i = 1/2 is a simple characteristic value of L, we may refer to [22, Thm. 1.40] to deduce that £ is the 
union of two subcontinua £ ± , where £ + consists of positive solutions near (1, 0), and £ + (and also £~) 
meets (1, 0) and either meets infinity in K x Ei or a point (/t + 1 /2, 0) with /i being a characteristic value 
of L different from 1/2. Consequently, £+ D (M + x E+) ^ 0, and we now show that £+ leaves K+ x E^" 
only at the bifurcation point (n, u) = (1,0). For, suppose the continuum £ + leaves R + x E^ at some point 
(n*,u*). Then there are (rij,Uj) eC + n (R + x E+) with 

(rij f Uj) — > (n*,u*) in RxEi . (2.22) 

In particular, writing rij = Aj + 1/2 we have itj = XjLuj + H(Xj, uf), j G N, and letting j — ► oo we 
obtain from Lemma [231 that u„ = A„Lu„ + -ff (A*, u») for A* := n* — 1/2. Hence 

9 a it* + A(n«, a)u» = , u*(0) = ^^^(u*) 

from which 

w* = n u » (•, 0)u*(0) , «*(0) = n*(3 Ut u*(0) 
by ( 12.191 ). Therefore, either u* (0) = and then = or (0) G E^~ \ {0} by ( 12.2b in which case > 
must be a characteristic value of Q Ut with a positive eigenvector. Thanks to ( 12.211 ) we derive 

u* =0 or u*(Q) G int(S+) . (2.23) 
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First suppose that = 0. Then (nj,Uj) — > (n*,0) in R+ x E+ , and we claim that A* = n* — 1/2 is a 
characteristic value of L. Indeed, putting vj := \\uj\\^Uj and taking into account ( 12.15b . the compactness 
of the operator L, and the fact that ||iT(A, m)||ei = °(IMIei) as MIei ^*0we may extract a subsequence 
of (vj) converging in Ex toward some v G \ {0} with v — \*Lv. Thus A* is indeed a characteristic 
value of L. As in ( 12.17b this implies that n* = A* + 1/2 is a characteristic value of Q with a positive 
eigenvector v(0) G E+, whence n* = 1 by assumption. Therefore, £ + leaves R + x Ef at the bifurcation 
point u») = (1, 0). Now suppose that 11,(0) G int(E+). Since £ + is connected and leaves R + x E^ at 
(n*, there is a sequence G £ + with Uj ^ E^ and (n^Uj) — > (n*, u«) in R x Ex. According 

to ( 12.2b we find m G N with u m (0) G E+. But (fi m , u m ) G £ + and thus 

d a u m + A(u m , a)u m = , u m (0) = n m i{u m ) , 

that is, u m = n Sm (• , 0)it m (0) G E+ due to ( 12.19b contradicting the choice of the sequence (Uj ) . Therefore, 
it* (0) G int(E+ ) is impossible. We have thus shown that £ + leaves R + x E^ only at the bifurcation point 
(n*,u*) = (1,0). 

It remains to prove that £ + does not meet a point [fi, + 1/2,0) with fi ^ 1/2 being a characteris- 
tic value of L. For, suppose £ + meets such a point. Then we find (rij,Uj) G £ + C R + x E^ with 
(nj,Uj) — > (/t + 1/2, 0) in R + x Ef . Exactly as above one shows that then /t = 1/2 which is ruled out by 
assumption. This proves the theorem. □ 

If the interior of the positive cone is empty, one can prove another result if we put some symmetry 
conditions on A and I. For the eigenvector B of Qq we only assume that B G Eq. 

Theorem 2.8. Suppose (T2~3Tl-(l2~9l), ( 12.19b , and ( 12.20b . Moreover, suppose that for each u G Ex, any positive 
eigenvalue ofQ u has geometric multiplicity 1 and possesses a positive eigenvector B u G Eq . In addition, 
fet A»(- u, •) = A*(u, •) and £*(— it, •) = £*(«,•) /or it G Ex. Ifr(Qo) — 1 is a simple eigenvalue 
°f Qo G K+(E q ) with an eigenvector B G Eq , ant/ if there is no other eigenvalue of Qq with positive 
eigenvector, then 

£+ := { (n, u) G £ ; tt(0) G } U { (n, -u) ; (n, u) G £ , u(0) , n > 0} 

ts an unbounded closed subset o/R x E^~ sttc/t f/zaf £ + \ {(1, 0)} consists of positive nontrivial solutions 
to ( 12.1b . where £ is the maximal continuum from Proposition \2.5\ 

Proof. As in the proof of Theorem l2.7l we have u G E^ for (n, it) G £ provided u(0) G Eq . If (n, it) G £ 
is such that tt(0) g" and n > 0, then u = H u (-, 0)tt(0) and u(0) = nQ n it(0). Hence n > is a 
characteristic value of Q« and by assumption there is a corresponding eigenvector B u G -Eq such that 
u(0) = r n B u for some r u < 0. Due to the symmetry conditions put on A* and it is easily seen that 
v := —it satisfies v = H v (-,0)v(0) and u(0) = nQ„i;(0) with v(0) = —r u B u G Eq~ . Consequently, 
(rt, i>) G R + x E+ solves (12.11 ). Hence £ + consists of nonnegative solutions only. To show that £ is 
unbounded assume to the contrary that £ meets a point (fi + 1/2, 0) with a characteristic value /t of L 
different from 1/2. Let (Aj + 1/2, Uj) G £ be such that (Xj,Uj) — > (//, 0) in R x Ex. Recall ( E^ and 
set := if itj(0) G .Eg" and := — Uj if 11^(0) ^ Eq. Then Uj G Ef and (Xj,v,j) — > (/t, 0) in 
R x Ex . Exactly as in the proof of Theorem l2.7l we deduce that /t + 1/2 is a characteristic value of Qo with 
an eigenvector in E+, whence ft + 1/2 = r(Qo) = 1 by assumption contradicting /t ^ 1 /2. Therefore, £ is 
unbounded according to Proposition ^ .5 1 and so is £ + . That the latter is closed in in R x Ex is a consequence 
of the continuity of L and H. This concludes the proof. □ 

Let us consider the set of possible parameter values in more detail. For, suppose the conditions of 
Theorem l2.71 let £ + denote the unbounded continuum in R + x E^ consisting of positive solution to ( 12.1b . 
and recall that Q u G K.+ (E q ) for u G Ex was defined by 



Q u w := £ o (n„(-,0)w) +e4u,U u (-,0)w) , weE,, u G Ex 
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Observe that for (n, u) 6 € + we have u = IT„(-, 0)w(0) with u(0) — nQ u u(0), whence nr(Q u ) > 1 for 
(n, ii) £ £ + . In addition, if 

for each m£E|,Q m e /C+(E ? ) has only r(Q u ) > as eigenvalue with positive eigenvector , (2.24) 

which holds e.g. if the Krein-Rutman theorem applies to Q u , then necessarily 

w(Q„) = l, (n,u)e£+. (2.25) 

This observations guarantees a more precise characterization of the spectrum 

a := {n; there is ueE| with (re,u) eC + \ {(1,0)}} 

as well as of the solution set 

T := {u; there is n G R with (n,u) eC + \ {(1,0)}} . 

The next proposition is in the spirit of [9| for the spatially homogeneous case A = in ( 11.31 ) and the 
easy proofs carry over to the present situation almost verbatim. We nevertheless include them here for the 
reader's ease. 

Proposition 2.9. Suppose the conditions of Theorem \2. 7\ and ( 12.241 . Setting 

Oi := inf a , a s := super , N t := inf r(Q u ) , N s := supr(Q u ) , 

we have: 

(i) ^ d C K+ andY C E+. 

(ii) < <7i < 1 < <7 S < oo. 

(iii) Ni — (^"cr s = oo, ami j/iVj > 0, then a s — 1/Ni < oo. 

(iv) iV, = oo iff (Ti — 0, and (/A^ < oo, then Oi — 1/N S > 0. 

(v) If r{Q u ) — > as |M|ei — * oo, f/ten cr s = oo. 

(vi) Ifr(Q u ) < £ for some £ > and every solution (n,u) G R + x E^ to J2.lt . f/ien er^ > l/£. In 
particular, ifr(Q u ) < 1 for every solution (n 7 u) G R + x E^~ to ( 12.11 ), f/ien er^ = 1 corresponding 
to supercritical bifurcation. 

Proof, (i) We note that, according to ( 12.19b . any solution (n, u) to ( 12. It satisfies it = n u (-,0)u(0) with 
w(0) = if n = 0, whence it = in this case. Hence cr. 

(ii) Since (1, 0) £ £ + this is immediate. 

(iii) Let Ni = 0. Then there is a sequence ((rij , Uj)) in £ + with Uj ^ and r(Q u - ) \ 0. From ( 12.251 ) 
we obtain rij /* oo, whence a s = oo. Conversely, let a s — oo. Then there is a sequence ((rij, Uj)) in 
£ + with rij / oo and ( 12.251 ) ensures r(Q Uj ) \ 0, whence iVj = 0. Let now N > 0. On the one hand, 
we necessarily have 1 = nr(Q u ) > nNi for all (n, u) G € + . Consequently, n < 1/Ni for n G a so that 
er s < 1/Ni. On the other hand, since 1 = nr(Q u ) < a s r(Q u ) for (n, u) G £ + we have r(Q u ) > l/cr s for 
it G T and thus N t > l/a s . 

(iv) The same as in (iii). 

(v) Theorem 12.71 implies that a s = oo or there are Uj G T with HmjHei — > oo. Then a s — oo or 
r{Qu) — > by assumption, hence er., = oo in both cases according to (iii). 

(vi) From (iv) we obtain <jj = 1/N S > l/£ if iV s < £. In particular, if £ = 1, then a; = 1 due to (ii). 

□ 
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Of course, particularly interesting is the case a s = oo. As in |9, Cor.3] one can easily impose conditions 
on b and A = A + \i guaranteeing r(Q u ) — > as ||u||ei - * 0, whence a s = oo by (v) of the above propo- 
sition. In particular, as in l28l Ex. 3. 3] one can put conditions on the same quantities leading to r(Q u ) < 1, 
that is, to supercritical bifurcation in view of (vi). 

As pointed out in Remark |2~6l the assumptions on the spectral radii and the properties of the eigenvalues 
of the operators Q u in Theorem l2.7l or Theorem l2.8l are not too restrictive but rather natural in applications 
due to the fact that these operators are compact and strongly positive in many cases. In general, we refer to 
ll28l Thm.3.1, Ex. 3.1, Ex. 3. 2, Ex. 3. 3] for other examples of diffusion operators A, birth moduli b, and death 
moduli p, satisfying the assumptions of Theorem |2.7| or Theorem |2.8| and provide here merely one example. 
Clearly, the set of applications is not restricted to following example. 



Example 2.10. Let Vl C M. , N > 1, be a bounded and smooth domain lying locally on one side of dfl. 
Let the boundary dQ be the distinct union of two sets Tq and T\ both of which are open and closed in dQ. 
For simplicity we assume a m G (0, oo) and consider a second order differential operator of the form 

A(u, a, x)w := — S7 X ■ (D(a, x)\7 x w) + g(u(a), V x n(a)) • V x w + h(u(a), V x u(a))w , (2.26) 

where 

D : J — > C 1 ^) is bounded and uniformly Holder 

_ (2.27) 
continuous with D(a, x) > , (a, x) £ J x f2 , 

and 

g e C 3 (R x R N , R N ) with g(0,0) = and h e C 3 (R x R N ,R) . (2.28) 

Let 

i/oeC 1 ^), u (x)>0, x GTi , (2.29) 
and let v denote the outward unit normal to Y\. Let 

w , on To , 

tj^ w : .-I I \ .1 hi - J. 1 . 



B{x)w 



J^w + v (x)w, onT 

Fix p, q £ (1, oo) with 



1 N 

- + — < 1 , (2.30) 
P 2q 

and let Eq := L q := L q (Q) be ordered by its positive cone of functions that are nonnegative almost 
everywhere. Observe that 

Ei ~ Wis := Wl B (Q) := {u e W 2 q ; Bu = 0} «-» L q = E 

and, up to equivalent norms, the interpolation spaces are subspaces of the Besov spaces B 2 ^ p := B 2 ^ p (il), 
that is, 

{ B 2 q %, 0<2£<l/g, 
E 6 := (L q , Wis) = B 2 q ls ■={ { w G B l% ; "lr = 0} , l/q < 2£ < 1 + l/q , 2^ 1 , 

fin <= R 2 « • Rn = f)\ 1 -\- ~\ / n <? 9.£ <? 9. 



{w G B 2 q % ;Bu = 0}, 1 + l/q < 2£ < 2 

? 2-2/j 



(see e.g. l|25l ). In particular, since 2(1 — 1/p) > N/q, we have E^ = B 2 2 g C(Q)forc; = 1 — 1/p 



awo? f/iMS int(_E+) 7^ 0. Assume then further that 

\i : R x J — > [0, 00) is uniformly Lipschitz continuous on bounded sets (2.31) 
and if 9(a) := /i(0, a) + h(0, 0), then 

9 : J — > (0, 00) ;s bounded and uniformly Holder continuous . (2.32) 
Finally, let b be such that 

b e C 3 (R, (0,oo)) . (2.33) 
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Consider 

A(u, a)w := A(u, a, -)w , w G E\ , ugEi, 
one/ define A(u, a) := A{u 1 a) + fi(u(a), a) for «£Ei and a G J. Then 

A a (a)w := A(0, a) to = -Vj, • (D(a, -)V x w) + 9{a)w 

is such that Ao : J — > C(W q g, L g ) is bounded and uniformly Holder continuous by (12.271 1 one/ (12.32b . 
Moreover, from [2, Sect.7,Thm. 11.1] it follows that for a G J fixed, — Ao(a) resolvent positive, gener- 
ates a contraction semigroup of negative type on each L r (Q), r G (1, oo), one/ is self-adjoint on L2{£1). 
Hence Ao generates a positive parabolic evolution operator with regularity subspace E\ — W q B by [4, 
II.Cor.4.4.2] one/ possesses maximal L p -regularity according to |4 III.Ex.4.7.3,III.Thm.4.10.10], whence 
(|23I >. Owing to d230b we may choose numbers s, s, and a such that 

- <s < s <l-a <~ + ^- - — . (2.34) 
p 2 2p 4q 

Tnen, as in Remark \2.1\ we have 

Ei := L p ( J, W£ B ) n J, L 9 ) — ( J, B ? 2 « B ) =: , j = 1, . . . , 4 , 

and analogously to 11261 Lem.2.7] we obtain that 

[u ^ g{u, V x u)] : W£(J, B£. B ) -> W^(J, B^ 1 ) 

is Lipschitz continuous. In particular, since s > 1/p and 2a — 1 > N/q we deduce for u G Ei fnaf 

[a i-> a(u(a), V x u(a))] : J -> C(Q) 

is Holder continuous. Clearly, the same holds true for h(u, V^u) ana! a/so 

[m : Wp(J, B 2 q a p . B ) -> Wp(J, B 2 q a p . B ) is Lipschitz continuous . (2.35) 

Similarly we obtain from (12. 3U and the embedding Ei C ?_w ( J, u < (being due to the interpo- 
lation inequality |4] I.Thm.2.1 1.1]) fnaf [a i— > /i(u(a),a)] : J — > C(f2) is Holder continuous. Gathering 
these information and invoking [4, II.Cor.4.4.2] we deduce that A(u, •) generates for each u G Ei a positive 
parabolic evolution operator Tl u (a, a), < a < a < a rn , on E — L q , whence d2.19l l. Obviously, ( 12.6l l 
holds for 

A*(u, a) := A(m, a) — Ao(a) . 
Afofe then that pointwise multiplication 

R 2a R 2a R 2(l-l/p) _j_ „ 

is continuous according to (12.341 anaf 13, Thm.4.11. 77z«s, the same theorem together with s > 1/p ensures 
that pointwise multiplication 

w;{j,b\%. b ) ■ w;{j,Bii. <B ) - ii(J,< ( ; ;e 1/p) ) 

is continuous since J — (0, a TO ) wifn a m < oo. From fnis we obtain £q G jC(Ei, £7i) ant/ ^* G C(F4, -E^), 
wnere 

£ (u):=b(Q) u(a)da and l*(u) := / [&(u(a)) - 6(0)]ti(a) da . 
Jo Jo 

/•Mrfnermore, 

< -6(0)]] W a WB 2« ;B) ||u|| w s (JiB 2« s) , 

whence d2.8l > anaf d2.9l > s/nce s > s. Defining 

b(u(a))IL u (a, 0) da , u G Ei , 

if is immediate that Q u G /C_|_(-B s ) awe to ||4] II.Lem.5.1.3] ana" fne compact embedding E\ c ^-> _E S since 
a m < oo. Moreover, Q u is strongly positive since II u (a, 0) is strongly positive on E^ for a G J \ {0} (see 
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Bill Thm. 13.6]j. In particular, Q u is irreducible and so, by 111 1 1 Thm. 12.3], r(Q u ) > is simple and the 
only eigenvalue of Q u with a positive eigenfunction since mt(E^) ^ as observed above. This ensures 
(12.20t and ( I2.2U . Therefore, ifb is normalized such that r(Qo) = 1, we conclude thanks to Theorem \2.7\ 

Proposition 2.11. Suppose ( 12.261 >-( |2~3"3T > and let r(Qo) = 1. Then the problem 

d a u + A(u, a, x)u + fj,(u(a), aju = , a£ J, a? G fi , 

/>oo 

u(Q,x)—n / b(u(a, x))u(a, x) da , x G f2 , 
Jo 

B(x)u(a, x) — , a>0, x £ dfl , 

admits an unbounded contiunuum £ + of positive nontrivial solutions (n, u) in 

M+ x (L+(J,^ B )nff;(J,L,)) . 

A noteworthy variant of the previous example is considering a functional (instead of a local) dependence 
of A, b, or p on u with respect to age. For details we refer to [28, Ex. 3. 2, Ex. 3. 3]. 



3. Positive Solutions via a Fixed Point Argument 

The aim of this section is to give sufficient conditions for the existence of nontrivial nonnegative solutions 
to the (parameter-independent) problem 

d a u + A(u, a) u = , a e J . 

(3.1) 



= / b(u,a)u(a)da , 
Jo 

in Eq without assuming a decomposition (I2.4b - d2.6b . Due to the quasilinear structure of the first equation we 
require some assumptions that can considerably be weaken if one restricts to linear problems. For 9 £ (0,1) 
we put X e := Li( J, Eg) and X@ := J, Eg). Let 

< a < j3 < 1 . (3.2) 

We suppose that, given any R > 0, there are p, u>, r\ > 0, a £ K, and k > 1 depending possibly on i? such 
that for $ a := B Xq (0, i?) n X+ we have 

A e C°' P ($ Q x J,£(Ei,E )) with 

[A(tv)] PfJ <»7, £r + A(«,-)cW(Bi,£b;K,«), (3.3) 
and A(it, •) is resolvent positive for each u e & a . 

Observe that d3 ,3b and [4, II.Cor.4.4] ensure that for each u G $ Q there is a unique positive parabolic 
evolution operator H u (a, er), < a < a < a m , on Eq corresponding to A(u, ■). The evolution operator 
satisfies according to |4 II.Lem.5.1.3] the estimates 

\\nu(a,a)\\ c{Es) + (a~ a)Z-V 2 \\n u (a,o-)\\ c{EoEe) < c e^ a ~^ , < a < a < a m , (3.4) 

for < C < £, < 1 ar, d some constants c = c (i?, £, C) > and ^ = v(R) G R (independent of u G <J> Q ). 
We assume that 

v < if a m = oo . (3.5) 

To control the dependence of the evolution operator n u on u G $ Q we require for each u G $ Q the existence 
of e = e(w) > and a measurable function g : (0, e) x J — > R+ (depending possibly on u) with 

lim / <7(r, a) da = , 

r-0+./o (3.6) 

max ||A(w,cr) - A(u,<j)\\ c{E E) < g(\\u - u\\ x +,a) , aeJ, \\u - u\\ x + < e . 

0<a<a a a 
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As for b appearing in ( 13. It we suppose that 

be E+($ Q xJ,F), 

(3.7) 

e va \\b(u,a) — b(u,a)\\F < g(\\u — u\\ x + : a) , aeJ, \\u — u\\ x + < e . 

Here F is assumed to be a Banach space ordered by a convex cone F + such that a (bilinear) multiplication 
m := [(/, e) i— > /e] is induced which is continuous considered as mappings 

m : F x E p -> E (j and m : F x E x -> E s for some S 6 (/?, 1] (3.8) 

and such that m(f, e) = fe G -Eg for / <E E+ and e 6 E^. Note that F = K is appropriate with 5 = 1. 



As a consequence of d3.41 >. (13 5b . (13.7b . and the compact embedding S<5 c ^ E^ we have 

Q(u)^ f m b(u,a)U u { ai O)daeC + (E p ,E s )nlC + (E p ), u £ $ Q . (3.9) 
Jo 

Solutions to (13.1b are, as noted in the introduction, fixed points of the map 

(u,B) » (II„(-, 0)B,Q(u)B) 
with B = u(0). Clearly, (13.2b - d3~8l are technical but not restrictive assumptions for applications (see 



Sect. 3]). However, since the main task is to single out nontrivial solutions we also have to impose structural 
and thus more restrictive assumptions in order to apply a fixed point theorem for conical shells (TJ. The 
assumptions read: 

there are r > and V £ E+ with tp £ (J rg + (l - Q{u)) , (3.10) 



\\v\\Xa <T 



where rg + (l - Q{u)) := {(1 - Q(u))B ; B e E^}, and 

+ 



there is t\ > such that r(Q(u)) < 1 for u e X{ with \\u\\ Xfj > n , (3.11) 
where r(Q(u)) denotes the spectral radius of the operator Q(u) G C(Ep). 

We comment in more detail on the structural requirements ( 13.10b . ( 13.1 lb after the proof of the following 
result, which is in the spirit of ll20l Thm.l]: 

Theorem 3.1. Suppose ( 13.2b , (13.3b , (I3.5b -( l3.8b , ( 13.10b , and (13.11b . Then (13.1b has at least one nontrivial 
nonnegative solution 

u e £i (J, Ex) n C l { J \ {0}, E ) n C(J, E s ) . 

Proof. We shall employ Q Thm.12.3] in proving the statement. Let X := X+ x E^ and X R := B x (0, R) 
with 

R ■= n(\\i\\c(Et,,E a ) + II&IU«(B v+ (0,t 1 )xJ,.F)) > , 

where i is the natural injection Ep E a . We put 

f(u,B) := (n B (.,0)B,Q(u)B) 

and first claim that / : X.r — > X is continuous and /(X.r) is relatively compact in X. Indeed, given 
(u, B), (u, B) £ Xh we note that 

/ ||n u (a, 0)B-U u (a, 0)B\\ Ea da < / ||n u (a, 0) - U u (a, 0)\\ C (E a ) da \\B\\ Ea 
Jo Jo 

+ f m ||n fi (a,0)|| £(££>) da||E-E|| Bci 
Jo 

<c{R) I g(\\u -u\\ x +, a) da + c(R)\\B - B\\ Ea , 



14 



CH. WALKER 



where we invoked S7~fy-SJM, and II.Lem.5.1.4]. Thus ILj(a,0)B -> U u (a,0)B in X a as (u,B) 
approaches (u, B) in by S3. 61 . Similarly we deduce Q(u)B — > Q(u)B in U^g as (m, B) — > (w, B) in 
Xfl, whence the continuity of /. Next, we use the characterization for compact sets in X a = L\(J, E a ) 
due to [14]. We may assume a m = oo. The previous argument entails 

sup ||n u (-,0)B||x Q < oo . 

(«,B)ex B 

Moreover, 



/ ||n u (a,0)B|| iSc , da < c(R) / e 1 " 1 da — >0 as 
Jn Jn 



N 



uniformly with respect to (u, B) G X^ by ( 13.4-b and ( 13.5b . Let ft, > 0. Then, from ( 13.4b , ( 13. 5b , and equation 
(II.5.3.8) in El we deduce 

p OO /"OC 

/ ||II (o + /i,0)B-n u (o,0)fl|U«,do < / ||n u (a + /i,0)-n u (a ! 0)|| £(E(3>B(i) da||B|U (3 
Jo Jo 

< c(R) hP- a 

and the right hand side tends to as h — > uniformly with respect to (u, B) 6 X# in view of ( 13.2b . 
Furthermore, since (13.4b . (13.5b ensure 

\\U u (a,0)B\\ El3 < c(R) , oe(0,oo), (u,B)€X fl , 

we obtain from the compact embedding Ep E a that IT u (a, 0)B belongs to a fixed compact subset of 
E a . Applying now [14, Thm.A.l] we derive the relative compactness of the set {n. u (-, 0)B ; (u, B) G X#} 
in X a . Next, observing that 



\\Q(u)B\ 



Es 



< 



\\b{u,a)\\ F \\n u {aM\c{E p , El )<ia\\B\\ Efi < c(R) 



for (u, B) G X# according to ( 13. 4b . (13.7b . and (13.8b we may use the compact embedding Ba CJ ^ B^ 
and also obtain the relative compactness of the set {Q(u)B ; (u, B) e Xjj} in B^. Therefore, /(X#) is 
relatively compact in X. It remains to check the crucial conditions (i) and (ii) from [1 1 Thm. 12.3]. For (i) 
suppose there exist A > 1 and (u, B) G X# for which 

\H\x ol + \\B\\E e = \\(u,B)\\ XR = R and f(u, B) — X(u, B) , 

that is, 

Xu(a) = II u (a,0)B , a G J , 
AB = Q{u)B . 

Since A > 1 we have B/0 (otherwise u = contradicting B > 0). From dl2i >. d3~4l . (|33I >. ( T3~9l l we 
deduce, on the one hand, that B G B^ is an eigenvector for Q(u) corresponding to the eigenvalue A > 1 
and, on the other hand, that u G X±. Invoking (13.11b we see that this is only possible if ||u||x 3 < T i- 
Consequently, recalling (13.8b and (13.12b we derive the contradiction 

|«k» + W b Wev = / \\u(a)\\ Ea da 



(3.12) 



R 



A 



&(«, a)II u (a, 0)B da 



< n INI 

< n||i|U(j5fl,B Q 



b(u, a)u(a) da 



r+ (0,n)xJ,F) 



< R. 



This ensures f{u, B) ^ X(u, B) for all A > 1 and all (u, B) G X# with || (it, B)||x B = B, whence (i) from 
0] Thm. 12. 3]. Finally, let ip be as in (13.10b with tq < R and assume there exists A > and (u, B) G X# 
with \\{u,B)\\ Xr =t and (u,B) — f(u,B) = X(0,ip). Thenu = ]!„(•, 0)B, hence u G Xi by ( H41 (|33T > 
with ||u||x a < tq, and B = Q(u)B + Xift. The latter implies -0 G rg , (1 — Q(u)) contradicting (13.10b . 
Thus (ii) from [1 , Thm. 12.3] is verified, too, and we conclude a fixed point (u, B) G X^ \ {(0, 0)} of the 
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map /, that is, a nontrivial positive solution to d3.lt . As for the additional regularity stated in the theorem 
we observe that necessarily u(a) — H u (a, 0)B for a € J with B = Q(u)B £ E^. It thus suffices to refer 
to the regularity theory of Chapter II in (4J- D 

Example 12. 101 or the examples in ||28l Sect. 3] apply with minor modifications to the situation of Theo- 
rem 13.11 We note that the special assumptions d3.5l l, ( 13.1 01 > , and (13. lit are also not too hard to verify in 
applications in view of the following remark. 

Remarks 3.2. (a) Let A(u, a) be of the form (11.31 > with /i being real-valued so that its evolution operator 
is given by 

U u (a,a)=e-S>^ dr U A(u , ) (a,a), 

where Um u a denotes the parabolic evolution operator corresponding to A(u, •). Then ( 13.5b holds provided 
that there is fig > such that lim fi(u, a) > /i uniformly with respect to u € $ Q and /io > s(— A(u, a)) 

a — >oo 

for u G $ Q and a G (0, oo) with s(—A(u, a)) being the spectral bound of the operator —A(u, a) consid- 
ered as a linear operator in Eg (see (4], Sect.I.l, Sect.II.5]). 

(b) Suppose (13.91 >. Then condition d3.1 lb is equivalent to assume that ker(A - Q{u)) nE| = {0}/or 
all A > 1 and u G X^~ with \\u\\xg > ti. 

Proof. This follows from (13.9b and the Krein-Rutman theorem which states that r(Q(u)) > is an eigen- 
value of Q(u) G 1C+(E[}) with a positive eigenvector. □ 

(c) Suppose (13. 9b . 77ie« condition Cl3 - 1 11 > ZioWi ||Q(iO||/;(.E e ) < 1 for some 9 G [0, (5] and all u G X^" 
wz'f/z [[tt||xja > r l- 

Proof. This is a consequence of (b) and (13.9b . □ 

f of) Suppose ( 13.9b . 77zen condition ( 13.10b is satisfied provided 

Q{u) - 1 G r+Cfy) /or u G X+ \ {0} tvftft ||u||x„ < r . 

Note that the latter condition corresponds in the non-diffusive case A = to assuming the scalar inequality 
Q{u) > lfor \u\ small as in l20l Thm.l] and [29, Thm.4.1]. 

Proof. Since in this case rg + (l — Q{u)) C — Et we may choose ip G Et \ {0} arbitrarily. □ 

(e) Suppose (13. 9b . Then condition d3.10b holds provided there is ip G ker(l — Q(u)) (~l Ej \ {0} such 

that Q{u) G K.+ (Ep) is irreducible for each u G X± \ {0} with \\u\\x a < r o fln <i <?•<?•> ^ e interior of Et 
is nonempty. 

Proof. If ip is as in the statement, then the Krein-Rutman theorem (e.g. see [jTT] Thm.12.3]) warrants that 
r(Q(u)) = 1 is a simple eigenvalue of Q(u), hence 

ker(l-Q(u))nrg(l-Q(u)) = {0}, 

from which we conclude ip g" rg , (1 — Q(u)) for each u. □ 

(f) If e .ff ects of small populations are negligible, then (13.101 ) holds. More precisely, let Q{u) — Q(0) for 
small u G and let Q(0) G K,+ {Ep) be irreducible, r(Q(0)) = 1, and mt(£7g ) 7^ 0. Then there is 
ip G Eg \ {0} vv/f/t "0 G ker(l — Q(0)) according to the Krein-Rutman theorem 111 II Thm.12.3], and ( 13.101 ) 
follows from (e). 
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(g) Suppose ( 13. 91 ). Then condition J 3 . 1 Ob is satisfied if there are To > and ip G int(E^) such that, 
given any u G X± \ {0} with ||u||xoi < T o> Q{ u ) !S irreducible and there exists a(u, •) G Li^ oc (0, a m ) 
with 

-A(u, a)ip = a(u, a)ip , aGJ, a«c/ / 6(u, a)e-^ a Q ( u < r ) dr da = 1 . 



o 



It is worthwhile to remark, however, that in this case a nontrivial solution to (13. It can be found also in the 
form u(a) = <j){a)ip, where the existence of a nonnegative nontrivial (j) follows from [20|. 



Proof. The assumptions imply n„ (a, O)^ = eti a M dr ^, aGJ, for u G X+ \ {0} with \\u\\ Xa < t , 
whence -0 G ker(l — Q(u)), and we may apply (e). □ 



(h) Clearly, Theorem \3.1\ applies to models involving several species, say with densities Uj, 1 < j < N, 
and u = (tii, ■ ■ ■ ,un)- If A and b in ( 13.11 ) have "diagonal form", that is, if each Uj satisfies ( 13.11 ) with 
A(u, a) and b(u, a) replaced by Aj(u, a) and bj(u, a), respectively, then it suffices to assume ( B.lOl l and 
( 13.1 11 1 for some component. More precisely, the assertion of Theorem 13.71 holds true provided there are 
j G {1, . . . , N} (at least one) and (13.10b is replaced by 

there are tq > 0, and ipj G E^with ipj g" r S+(l — Qj( u )) 

ue(x+r 

\\ u \\(X a ) N <T ° 

while d3.11| ) is replaced by 

there is n > such that r(Qj(u)) < 1 for u G (Xj^) with || , u||(x s ) JV > r i i 

Qj(u) := / 6j(m, a)n uj (a, 0) da 



o 

w/f/i n uj denoting the parabolic evolution operator corresponding to Aj(u, •). 

Proof. Looking for solutions u — (u±, . . . ,un), where only the j-components are non-vanishing, this 
follows by an obvious modification of the proof of Theorem |3.1| □ 
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